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CHAPTER  I 


INTRODUCTION 

General  Statement 

The  Raton  Basin  of  New  Mexico  and  Colorado  (35°-38°N;  104°-105.5°W) 
is  a  structural  depression  which  is  defined  by  the  Sangre  de  Cristo  Up¬ 
lift  on  the  West,  the  Wet  Mountains  Uplift  to  the  north,  the  Apishipa 
Arch  on  the  northeast,  and  the  Sierra  Grande  Arch  on  the  southeast 
(figure  1) .  Today  the  Raton  Basin  is  outlined  by  topographic  relief 
caused  by  the  caprock  "effect"  of  the  Trinidad  Sandstone.  The  importance 
of  the  outlining  "effect"  of  the  Trinidad  Sandstone  is  illustrated  by 
both  Baltz  (1965)  and  Brill  (1952)  who  draw  the  eastern  boundary  of  the 
Raton  Basin  at  the  outcrop  of  Trinidad  Sandstone,  rather  than  along  the 
axes  of  the  Apishipa  and  Sierra  Grande  Arches. 

Extending  175  miles  long  and  a  maximum  of  65  miles  wide,  the  Raton 
Basin  is  divided  by  the  Cimarron  Arch  into  the  northern  Raton  Basin  and 
the  Las  Vegas  Basin  on  the  south.  The  deepest  depression  within  the 
Basin  is  located  near  the  Sangre  de  Cristo  Uplift.  Rocks  on  the  west 
limb  of  the  Basin  dip  vertically,  whereas  rocks  on  the  east  limb  dip  one 
to  five  degrees  (Wanek,  1963).  In  this  thesis,  the  entire  Raton  Basin 
will  be  referred  to  as  the  Basin,  while  the  northern  portion  will  be 
called  the  Raton  Basin,  and  the  southern  portion  referred  to  as  the  Las 
Vegas  Basin. 
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Figure  1.  Geologic  structures  defining  the  Raton  Basin  of  New 
Mexico  and  Colorado. 


Purpose  of  the  Investigation 

Knowledge  of  the  Basin  is  incomplete,  and  correlation  of  geophysical 
and  geologic  data  has  not  been  previously  published.  Therefore,  in  order 
to  fill  this  void  in  the  literature  this  thesis  will  undertake  the  fol¬ 
lowing  objectives:  1)  Obtain  and  interpret  both  Free  Air  and  Bouguer 
Gravity  maps;  2)  Obtain  and  interpret  a  vertical  component  magnetic  map; 
3)  Obtain  vertical,  short-period  seismic  records  for  the  Capulin  area  in 
an  effort  to  determine  the  natural  seismic  activity;  4)  Compare 
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geophysical  data  with  the  geologic  structure,  and  5)  Determine  whether 
economic  potential  in  the  Basin  can  be  inferred  from  geophysical  data. 

The  study  being  reported  here  includes  gravity,  magnetic,  and  seismic 
surveys  of  portions  of  the  Basin,  including  the  Sierra  Grande  Arch  and 
adjacent  High  Plains  (figure  2). 

General  Geology  of  the  Basin 

The  first  of  the  two  orogenies  that  have  affected  the  Basin  occurred 
during  Pennsylvanian  and  Permian  time  (Baltz,  1965),  and  was  respon¬ 
sible  for  detritus  filling  the  Central  Colorado  Basin  (located  in  the 
northern  portion  of  the  Raton  Basin).  The  detrital  material  was  derived 
from  both  west  (San  Luis  Uplift)  and  east  (the  ancestral  Wet  Mountains 
Uplift)  of  the  Basin.  Contemporaneously,  the  Rowe-Mora  Basin,  which  was 
in  the  general  vicinity  of  the  present  Las  Vegas  Basin,  also  formed. 

This  Basin  was  defined  by  the  San  Luis  Uplift  to  the  west,  the  Sierra 
Grande  Uplift  to  the  east,  the  Cimarron  Arch  on  the  north,  and  was  linked 
by  a  low  saddle  to  the  Tucumcari  Basin  to  the  south.  The  Colorado  and 
Rowe-Mora  Basins  filled  with  sediments  which  were  derived  from  the 
ancestral  Rockies  and  deposited  in  a  marine  transgressive  environment  of 
Pennsylvanian  age.  The  resulting  strata  show  a  laterally  complex  compo¬ 
sition  reflecting  the  mixed  terrestial-marine  origin.  These  strata  are 
thickest  in  the  regime  of  the  Sangre  de  Cristo  Mountains.  Large  accumu¬ 
lations  of  sandstone  and  limestone  occur  in  the  southern  portion  of  the 
Rowe-Mora  Basin,  but  shale  deposits,  which  are  dominated  by  thick  grey 
to  black  strata,  predominate  to  the  north.  In  fact,  this  shale  which 
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characterizes  the  northern  part  of  the  ancient  Central  Colorado  Basin  is 
a  minimum  of  9,500  feet  thick  (Brill,  1952).  On  the  margins  of  the 
ancient  basins,  shale  deposits  thin  towards  the  San  Luis,  Sierra  Grande, 
and  Wet  Mountains  Uplifts. 

The  depositional  history  of  the  strata  in  the  Basin  is  complex  as 
is  the  correlation  between  contemporaneous  formations  in  the  Las  Vegas 
and  Raton  Basins  (figure  3) .  Strata  which  were  deposited  in  the  Permian 
era  are  dominated  by  marine  shales  and  clastic  rocks.  Triassic  strata 
are  composed  primarily  of  terrestial  clastic  rocks  but  are,  however, 
sparsely  represented  in  the  area.  Today  the  Las  Vegas  Basin  area  is 
covered  partially  by  Quaternary  volcanic  rocks.  Discussion  of  the  depo¬ 
sitional  and  orogenic  history  is  limited  to  the  area  of  the  present  Basin 
and  is  based  on  rock  outcrops  and  the  known  geologic  structure  (plate  1) . 

The  oldest  rocks  identified  in  the  Basin  are  marine  sandstones  and 
limestones  of  the  Devonian  Espiritu  Santo  Formation  (Baltz,  1965)  in  the 
Las  Vegas  Basin  and  the  Devonian  Chaffee  Formation  (Johnson,  1969)  in  the 
Raton  Basin.  Limestone  breccia  of  the  Mississippian  Tererro  Formation 
rests  unconformably  on  the  Espiritu  Santo  Formation  in  the  Las  Vegas 
Basin.  The  Tererro  Formation  is  missing  in  the  Raton  Basin. 

The  Pennsylvanian  Magdalena  Group  is  composed  of  sandstone  and  grey 
shale  of  the  Sandia  Formation  and  limestone  of  the  overlying  Madera  For¬ 
mation.  The  Sandia  Formation  in  the  Las  Vegas  Basin,  which  consists  of 
over  1,000  feet  of  sedimentary  rock,  is  correlative  with  the  Deer  Creek 
Formation  and  the  lower  part  of  the  Mlnturn  Formation  of  the  Raton  Basin. 
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Figure  3.  Columnar  section  of  the  Las  Vegas  and  Raton  Basins.  Formation 
thicknesses  are  from  Baltz  (1962);  symbols  used  in  the  columns 
are  from  Compton  (1962) . 
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The  Madera  limestone  ranges  from  1,000  to  3,000  feet  in  thickness  in  the 
Las  Vegas  Basin;  it  grades  laterally  northward  into  the  Mintum  Forma¬ 
tion,  which  has  a  thickness  of  5,500  feet.  These  units  thin  to  the 
west  and  east  of  the  Basin  in  the  general  direction  of  the  San  Luis  and 
Sierra  Grande  Uplifts. 

The  Upper  Pennsylvanian/Permian  Sangre  de  Cristo  Formation  is  ap¬ 
proximately  3,500  feet  thick  in  the  Las  Vegas  Basin  and  greater  than 
9,500  feet  thick  in  the  Raton  Basin  (figure  3). 

At  the  end  of  the  Permian  time  a  sea  filled  the  Rowe-Mora  Basin, 
and  deposition  of  siltstone  and  shale  of  the  Yeso  Formation  began. 

Located  above  the  Yeso  Formation  is  the  Glorieta  Sandstone  and  San 
Andreas  Limestone,  respectively.  Overlying  the  San  Andreas  Limestone 
is  the  Bernal  Formation  which  is  composed  of  limestone,  siltstone,  shale, 
and  sandstone.  These  three  formations  have  a  combined,  maximum  thickness 
of  1,075  feet  (figure  3). 

The  Late  Trlassic  Dockum  Group  lies  conformably  on  Permian  strata. 

In  the  Las  Vegas  Basin  the  Group  is  composed  of  both  the  Santa  Rosa 
Sandstone  and  the  overlying  Chinle  Formation  shales.  North  of  the  Cimar¬ 
ron  Arch  the  Dockum  Group  is  correlative  to  the  Johnson  Gap  Formation, 
which  is  composed  of  limestone  conglomerate,  limestone,  siltstone,  shale, 
and  sandstone  (Baltz,  1965).  The  Dockum  Group  is  approximately  1,100 
feet  thick  in  the  Las  Vegas  Basin  and  less  than  700  feet  thick  in  the 
Raton  Basin. 

By  Early  Cretaceous  time  the  region  was  part  of  the  Rocky  Mountain 
geosyncline.  The  Laramide  Orogeny,  the  second  orogeny  to  affect  the 
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area,  continued  from  Cretaceous  time  through  Oligocene  time  and  produced 
the  present  Las  Vegas  and  Raton  Basins  and  other  associated  structures. 

At  the  beginning  of  the  orogeny,  the  San  Luis  Uplift  was  still  a  positive 
area.  By  Eocene  time  the  present  Raton  Basin  was  defined  by  the  marked 
uplift  of  the  Sangre  de  Crlsto  region  and  the  gentle  uplift  of  the 
Apishipa-Sierra  Grande  Arches.  The  Sangre  de  Crlsto  Uplift  merged  with 
the  San  Luis  Uplift  prior  to  the  formation  of  the  Rio  Grande  rift.  When 
the  Rio  Grande  rift  formed  in  Miocene/Pliocene  time,  the  San  Luis  posi¬ 
tive  area  was  downfaulted  into  the  rift  graben. 

The  Early  Cretaceous  Purgatolre  Formation  is  composed  of  both  con¬ 
glomeratic  sandstone  and  sandstone.  The  Dakota  Sandstone  overlies  the 
Purgatolre  Formation  and  is  composed  of  interbedded,  buff  sandstone  and 
shale.  Together,  these  two  formations  are  up  to  650  feet  thick,  and 
today  they  form  the  primary  groundwater  aquifer  in  the  Basin.  Estimates 
place  the  volume  of  water  stored  in  the  Dakota  Sandstone  at  55  million 
acre-feet  in  New  Mexico  (Griggs,  1948);  twice  the  above  volume  is  esti¬ 
mated  for  the  Baain. 

The  Graneros  Shale,  Greenhorn  Limestone,  and  Carlile  Shale  are  Late 
Cretaceous  in  age  and  vary  in  thickness  from  385  feet  to  700  feet  (Baltz, 
1965) .  The  Upper  Cretaceous  Niobrara  Formation  is  composed  of  less  than 
100  feet  of  shale  and  sandstone  and  rests  conformably  on  the  Carlile 
Shale. 

Located  on  top  of  the  Niobrara  Formation  is  the  Upper  Cretaceous 
Pierre  Shale,  which  is  up  to  2,300  feet  thick  in  the  Raton  Basin  but 
is  absent  in  the  Las  Vegas  Basin.  The  formation  is  composed  of  shale. 


thin  beds  of  limestone,  and  sandstone.  Above  the  Pierre  Shale  are  the 
Upper  Cretaceous  Trinidad  Sandstone  and  Vermejo  Formation  which  are  un¬ 
differentiated  on  some  geological  maps  (Bachman  and  Dane,  1962).  The 
Trinidad  Sandstone  is  composed  of  arkosic  sandstone  with  thin  interbed- 
ded  shale;  it  intertongues  with  the  shale,  coal,  and  arkosic  sandstone 
of  the  Vermejo  Formation,  The  thickness  of  the  two  formations  varies 
from  250  feet  to  850  feet.  The  Vermejo  Formation  is  the  main  source  of 
coal  in  the  Basin,  producing  high  volatile  C  bituminous  coal  of  coking 
quality  (Wanek,  1963).  In  the  areas  of  extensive  intrusive  rocks,  much 
of  the  coal  was  destroyed  by  contact  metamorphism  (Jurle  and  Gerhard, 
1969). 

The  Raton  Formation  is  Late  Cretaceous  to  Early  Tertiary  in  age  and 
overlies  the  Vermejo  Formation.  It  has  a  maximum  thickness  of  1,700 
feet  and  is  composed  of  arkosic  sandstone,  shale,  and  coal.  The  Poison 
Canyon  Formation  of  Paleocene  age  lies  unconformably  on  the  Raton  Forma¬ 
tion.  It  consists  of  as  much  as  2,500  feet  of  arkosic  sandstone,  con¬ 
glomerate,  and  thin  shale.  The  youngest  significant  formation  in  the 
Basin  is  the  Cuchara  Formation, which  is  exposed  around  Spanish  Peaks. 

It  is  Eocene  in  age  and  is  composed  of  conglomeratic  sandstone  and  inter- 
bedded  shale. 

During  the  latter  stages  of  the  Laramide  Orogeny  numerous  intrusive 
rocks  penetrated  the  sedimentary  rocks.  Dike  swarms,  sills,  stocks  and 
laccoliths  dating  from  Eocene/Oligocene  time  are  found  throughout  the 
Basin.  However,  the  greatest  concentration  is  found  in  the  Spanish  Peaks 
intrusive  area.  Many  small  anticlines  are  attributed  to  the  intrusive 
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upvelling.  Emplacement  of  many  of  the  dikes,  which  are  silicic  to  ultra 
basic  In  composition,  appears  to  have  been  fault  controlled  (Wanek, 

1963). 

Volcanic  activity  was  concentrated  near  the  Sierra  Grande  Arch. 
Baldwin  and  Huehlberger  (1959)  identified  three  sequences  of  volcanic 
flows,  dating  from  Miocene  to  Recent,  that  they  called  Raton,  Clayton, 
and  Capulin,  respectively.  Both  the  Raton  and  Capulin  sequences  are  com¬ 
posed  of  basalts,  and  the  Clayton  sequence  has  a  varied  composition  which 
ranges  from  extremely  alkalic  to  subs  Hide.  Geochemical  analysis  of 
the  volcanic  rock  suggests  the  magma  originated  in  the  upper  mantle 
(Jones  et  al . ,  1976).  Stormer  (1972)  suggests  eastward  migration  of 
the  vulcanism  across  the  Rio  Grande  rift,  from  the  San  Juan  volcanic 
field  to  the  Sierra  Grande  volcanic  field,  occurred.  Sanford  et  al. , 
(1981)  also  suggest  that  the  Sierra  Grande  vulcanism  might  be  an  eastern 
extension  of  the  Jernez  Lineation  of  western  New  Mexico.  Further, 

Hamilton  and  Pakiser  (1965)  indicate  that  the  Rio  Grande  rift  is  limited 
to  the  upper  crust.  Thus,  the  mantle  may  be  the  magma  source  in  vulcan¬ 
ism  on  both  sides  of  the  rift. 

It  has  been  suggested  (Edwards  et  al.,  1978)  that  anomalously  high 
heat  flows  may  be  caused  by  igneous  activity  within  the  Basin.  Typical 

2 

heat  flow  values  for  the  Great  Plains  average  1.5  HFU  (1  HFU  *  lycal/cm 
sec),  but  values  as  high  as  4.7  HFU  are  found  in  the  Basin.  Suppe  et 
al. ,  (1975)  proposed  a  hot  spot  trace  from  eastern  Arizona  to  a  location 
near  Raton,  New  Mexico.  Reiter  et  al.,  (1979)  theorized  that  a  thermal 
anomaly  formerly  existed  east  of  the  Southern  Rocky  Mountains  and  is  now 
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centered  near  the  Sierra  Grande  Arch.  They  noted  that  radioactivity  of 
the  sedimentary  rocks  within  the  Basin  can  account  for  only  ten  percent 
of  the  observed  increased  heat  flow;  the  remainder  is  unexplained. 

In  summary,  the  total  thickness  of  sedimentary  fill  in  the  Basin  is 
about  17,000  feet.  The  Basin  sediments  date  from  Middle  Paleozoic  time 
and  have  been  affected  by  two  orogenies.  These  sedimentary  rocks  have 
economic  importance,  and  mining  activities  in  the  Raton  Basin  have  pro¬ 
duced  coal  and  graphite  in  commercial  quantities .  Cretaceous  rocks  are 
also  known  to  have  both  oil  and  gas  in  small  quantities.  Currently, 
there  is  one  commercial  oil  well  in  Huerfano  County,  Colorado,  at  the 
northern  end  of  the  Basin,  and  a  gas  field  near  the  Basin  in  eastern  Las 
Animas  County,  Colorado.  A  basement  structure  map  will  be  developed  in 
Chapter  III  using  the  densities  of  the  sediments  described  above,  and 
the  gravity  data  that  will  be  introduced  in  Chapter  II. 


CHAPTER  II 


DATA  COLLECTION  AND  REDUCTION 

General  Statement 

A  portion  of  the  gravity  data  used  in  this  study  were  obtained  from 
the  National  Oceanic  and  Atmospheric  Administration  (NOAA) ,  Washington, 
D.C.  The  gravity  and  magnetic  measurements  in  the  vicinity  of  Capulin 
National  Monument,  New  Mexico  and  the  magnetic  measurements  in  the  Raton 
Basin  were  made  by  the  author  during  1981.  Each  gravity  station  in  the 
vicinity  of  Capulin  National  Monument  was  also  a  magnetic  station.  In 
order  to  provide  a  detailed  analysis  of  the  Basin,  additional  magnetic 
stations  were  spaced  midway  between  the  gravity  stations.  NOAA  station 
density  is  0.06341  station  per  square  mile;  Capulin  gravity  station  den¬ 
sity  is  1.333  stations  per  square  mile;  and  Raton  Basin  magnetic  station 
density  is  0.0650  station  per  square  mile. 

Elevation  Control 

U.S.  Geological  Survey  fifteen  minute  quadrangle  maps  were  used  to 
determine  elevations  for  stations  near  Capulin.  Such  elevation  data 
were  supplemented  with  a  preliminary  Bureau  of  Land  Management  map  of 
Capulin  National  Monument.  Spot  elevations  on  the  maps  at  road  inter¬ 
sections  and  benchmarks  were  used  wherever  possible.  Elevation  accuracy 
is  within  five  feet,  corresponding  to  +0.35  mgal  accuracy  in  gravity 
measurement.  Elevation  effects  were  insignificant  in  the  reduction  of 
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the  magnetic  data. 

Magnetic  Data  Measurement  and  Reduction 
The  Capulin  magnetic  data  were  collected  using  an  E.J.  Sharpe  PMF-3 
flux  gate  magnetometer  (SO  and  100  gammas  per  division) .  Station  number 
one  (figure  14)  was  arbitrarily  chosen  by  the  author  as  the  datum  for 
the  survey.  Station  reoccupation  resulted  in  a  minimum  of  three  datum 
measurements  each  day  and  a  minimum  of  two  cross  correlations  on  stations 
occupied  on  previous  days.  Diurnal  correction  graphs  were  computed 
daily,  and  were  used  to  adjust  all  measurements. 

The  Raton  Basin  magnetic  data  were  collected  using  an  Askania  verti¬ 
cal  Torsion  magnetometer  (2.433  gammas/degree-scale  division).  Field 
procedures  for  the  Raton  magnetic  survey  were  the  same  as  those  employed 
in  the  Capulin  magnetic  survey.  Normal  field  corrections  were  computed 
from  36°30'N,  105°15'W  using  the  Regional  Vertical  Intensity  Map  1965. 
Latitude  correction  was  909.09  gammas  per  degree  latitude  (12.98  gammas 
per  mile);  longitude  correction  was  434.80  gammas  per  degree  longitude 
(7.9  gammas  per  mile). 

Gravity  Data  Measurement  and  Reduction 
Datum  for  gravity  data  acquired  from  N0AA  is  sea  level.  The  density 
used  in  Bouguer  anomaly  calculations  was  2.67  g/cc.  Gravity  data 
acquired  near  Capulin  National  Monument  were  collected  using  a  North 
American  gravimeter  (0.107  mgal  per  division).  Station  number  one  was 
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selected  to  coincide  with  a  NOAA  station  in  the  survey  area  and  was  used 
as  the  datum  for  the  Capulin  survey.  This  procedure  allowed  comparison 
of  the  Capulin  data  to  the  NOAA  data.  Latitude  corrections  were  made 
from  36°48'18"N  at  a  rate  of  1.25  mgal  per  mile. 

Seismic  Data 

A  short  period  vertical  Geotech  18300  seismometer  with  a  free  period 
of  one  second  was  buried  approximately  two  feet  in  depth  at  the  Visitor's 
Center,  Capulin  National  Monument,  from  16  July  1981,  to  23  September 
1981.  During  this  period,  51  daily  seismograms  containing  1,122  hours  of 
data  were  recorded  by  a  Sprengnether  VR-50  recorder-amplifier.  Damping 
was  increased  after  14  August  1981  from  normal  internal  resistance  to 
2,200  ohms  resistance  across  the  input  terminals  of  the  amplifier.  Fil¬ 
ters  were  not  used. 

Each  seismogram  contains  twenty-four  hours  of  data  with  the  minute 
marks  spaced  sixty  millimeters  apart.  The  ink  trace  was  approximately 
1.0  mm  thick  and  makes  it  impossible  to  resolve  frequencies  higher  than 
one  Hertz.  Estimated  amplification  from  artificial  inputs  is  between 
25,000  to  50,000. 


Computer  Processing  Techniques 

Modified  computer  programs  from  Davis  (1973)  and  new  programs  writ¬ 
ten  by  the  author  for  this  study  were  employed  to  process  both  magnetic 
and  gravity  data  on  an  IBM  360  computer  housed  at  Texas  Tech  University. 
The  machine  processing  helped  to  verify  manual  calculations  and  also  to 
derive  mathematical  models  that  would  have  been  impossible  to  obtain  by 
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other  means. 

All  data  were  systematically  processed  in  the  same  manner.  First, 
the  corrected  data  were  filtered  by  a  uniform  matrix  of  either  1,600  or 
2,400  elements.  The  matrix  value  of  each  element  was  obtained  by  calcu¬ 
lating  the  algebraic  average  of  the  six  nearest  data  values.  Development 
of  the  matrix  by  interpolation  over  six  data  values  for  each  element  re¬ 
sults  in  the  minimum  and  the  maximum  values  in  the  completed  matrix  to 
be  less  than  the  minimum  and  the  maximum  values  in  the  input  data.  Less 
than  six  data  values  in  the  interpolation  increased  the  minimum-maximum 
range  in  the  matrix;  however,  the  areal  extent  of  the  maximum  and  the 
minimum  features  would  have  appeared  exaggerated  in  the  final  map.  This 
causes  mountain  peaks  to  appear  as  elevated  buttes  or  mesas.  The  use  of 
more  than  s  i  data  values  in  interpolation  resulted  in  excessively  reduc¬ 
ing  the  minimum-maximum  range  in  the  matrix. 

Statistical  Analysis  System's  SAS/GRAPH  (Reinhardt,  1981)  was  used 
to  produce  the  computer  maps  from  the  interpolated  matrix.  These  maps 
were  contoured  on  a  strict  mathematical  ratio  by  the  computer;  therefore, 
the  technique  can  be  used  advantageously  to  compare  different  data  sets 
(gravity  and  magnetic)  produced  at  the  same  scale  and  in  the  same  geo¬ 
graphic  area.  In  order  to  compare  the  different  mapping  techniques, 
both  computer  plotted  and  hand  drawn  maps  using  the  same  data  sets  have 
been  included  in  this  report. 

The  computer  also  was  used  to  compare  the  computer  contoured  map  of 
the  data  contained  in  matrix  form  to  a  surface  defined  by  a  polynomial 
of  order  n.  A  polynomial  of  order  0  represents  a  dipping  plane,  and 
polynomials  of  very  high  order  very  accurately  represent  the  data  in 


matrix  form.  A  second  order  polynomial  was  found  to  accurately  represent 
the  steep  gravity  dip  under  the  Sangre  de  Cristo  Mountains  and  the  shal¬ 
low  gravity  dip  of  the  Great  Plains.  The  residual  gravity  (Smith,  1970; 
Dobrin,  1976),  or  difference,  of  the  second  order  polynomial  from  the 
matrix  gravity  data  was  used  to  construct  the  basement  structure  map 


described  in  Chapter  III. 


CHAPTER  III 


PRESENTATION  AND  INTERPRETATION  OF  THE 
GRAVITY,  MAGNETIC,  AND  SEISMIC  DATA 

Free  Air  Gravity  Data 

The  NOAA  free  air  gravity  data  have  been  tabulated  and  included  in 
Appendix  A.  The  gravity  data  acquired  in  the  Capulin  area  as  part  of 
this  study,  which  will  be  discussed  as  an  addition  to  the  NOAA  data,  have 
been  tabulated  and  included  in  Appendix  8. 

The  NOAA  free  air  gravity  anomaly  data  are  shown,  for  comparison,  in 
hand  drawn  (figure  4)  and  computer  plotted  (figure  5)  formats.  The  alge¬ 
braic  averages  for  the  station  elevation,  free  air  anomaly,  and  Bouguer 
anomaly  in  each  one  degree  block  in  the  study  area  are  shown  in  figure 
6.  The  free  air  anomalies  for  the  northern  two  blocks  of  the  study  area 
are  nearly  zero  and  suggest  isostatic  equilibrium  for  this  area.  How¬ 
ever,  the  free  air  anomalies  are  significantly  positive  in  the  southern 
two  blocks  of  the  study  area  and  could  indicate  lack  of  local  isostatic 
adjustment  in  these  areas.  The  question  of  regional  isostatic  equilibri¬ 
um  is,  therefore,  raised. 

Perfect  isostatic  adjustment  assumes  that  any  near  surface  positive 
mass  has  an  equal  amount  of  negative  mass  below  the  compensation  level 
within  the  Earth,  and  in  the  simple  case,  such  a  mass  distribution  would 
cause  near  zero  free  air  gravity  anomalies.  However,  in  elevated  ter¬ 
rain,  as  in  the  study  area,  the  gravitational  effect  upon  a  surface 
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Figure  4.  Regional  free  air  gravity  map  using  NOAA  gravity  data. 

Datum  is  sea  level  and  gravity  stations  are  shown  in 
Figure  10.  Contoured  by  hand. 
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6.  Graphic  representation  of  the  algebraic  averages  of  the 
NOAA  gravity  and  elevation  data  divided  into  one  degree 
latitude  and  longitude  blocks.  E  ■  elevation;  FA  “  free 
air  anomaly;  BA  -  Bouguer  anomaly;  FC  -  Faye  correction; 
EE  -  equivalent  elevation;  CFA  -  corrected  free  air 
anomaly  (FA  +  FC) . 
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measurement,  of  the  near  surface  rock  (the  positive  mass)  is  larger  than 
the  gravitational  effect  of  the  negative  mass  at  the  compensation  level 
simply  because  the  near  surface  rock  is  nearer  the  measured  station. 

This  means  that  a  positive  free  air  anomaly  is  in  part  caused  by  depth 
to  compensation  and  does  not  necessarily  indicate  crustal  loading.  It  is 
possible  that  the  total  observed  free  air  anomaly  can  be  explained  in 
terms  of  the  depth  of  isostatic  compensation.  Therefore,  if  the  free 
air  anomaly  is  to  be  used  in  an  estimation  of  isostatic  adjustment,  a 
correction  must  be  applied. 

The  correction  applied  to  the  observed  free  air  anomaly  for  the  ef¬ 
fect  of  depth  to  compensation  is  generally  called  the  Faye  correction, 
after  the  man  who  first  discussed  the  effect  in  the  Nineteenth  Century. 
One  method  of  calculating  the  Faye  correction  uses  the  Bouguer  anomaly 
to  determine  the  elevation  of  the  surface  of  a  slab  which  is  in  isostatic 
balance  (Woollard,  1962).  This  is  done  by  assuming  if  isostatic  adjust¬ 
ment  is  complete,  the  Bouguer  anomaly  is  related  to  the  surface  elevation 
of  the  crustal  load  causing  the  adjustment.  The  elevation  of  this  sur¬ 
face  slab  is  thought  of  as  the  equivalent  elevation  (h')  derived  from  the 
Bouguer  anomaly  (h'  ■  BA/2Tryp;  where  y  is  the  gravitational  constant  and 
p  is  the  assumed  density) .  The  equivalent  elevation  is  subtracted  from 
the  station  elevation  (h)  and  the  difference  (Ah)  defines  the  slab  thick¬ 
ness  used  in  the  Faye  correction  (F.C.):  F.C.  *  2TrypAh.  The  Faye  cor¬ 
rection  is  subtracted  from  the  free  air  anomaly  to  produce  a  corrected 
free  air  anomaly  value,  which  is  actually  the  isostatic  anomaly  since 
this  computation  has  removed  the  effect  of  depth  of  compensation  upon 
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the  free  air  anomaly.  As  a  result  of  applying  the  Faye  correction  to 
the  free  air  anomalies  in  the  study  area  (figure  6)  the  corrected  free 
air  anomalies  clearly  suggest  regional  isostatlc  equilibrium. 

Although  regional  isostatic  equilibrium  is  established,  the  lo¬ 
cally  positive  average  free  air  anomalies  of  the  southern  two  blocks 
(figure  6)  are  indications  of  a  lack  of  local  compensation.  The  Quater¬ 
nary  volcanic  rock  covering  most  of  the  southern  two  blocks  of  the  study 
area  is  missing  in  the  northern  portion  of  the  study  area  and  is  assumed 
to  be  the  cause  of  the  large  positive  free  air  anomaly  in  the  southern 
two  blocks.  If  this  is  so,  the  free  air  anomaly  can  be  used  to  estimate 
the  rock  thickness.  If  a  density  of  2.67  g/cc  is  assumed  for  the  volcan¬ 
ic  rock,  and  the  formula  for  the  gravitational  effect  of  a  slab  is  solved, 
a  thickness  of  180  m  is  indicated  for  the  Quaternary  volcanic  rocks.  As 
a  confirmation  of  the  rock  thickness  calculated  from  gravity  data,  a 
water  well  at  Capulin  National  Monument  penetrated  196  m  of  volcanic 
rock.  The  real  thickness  is  probably  somewhat  irregular,  and  the  value 
computed  from  the  gravity  data  appears  to  be  a  good  estimate. 

Comparison  of  the  free  air  map  and  topographic  map  (figure  7)  shows 
that  positive  topographic  features  are  also  areas  of  positive  free  air 
gravity  anomalies  throughout  the  Basin  and  Great  Plains.  This  relation¬ 
ship  is  to  be  expected  (Woollard,  1959),  and  the  larger  topographic  fea¬ 
tures  are  associated  with  larger  local  positive  free  air  anomalies. 

The  area  of  the  Capulin  free  air  gravity  map  (figure  8)  is  covered 
by  volcanic  rocks  of  three  sequences  -  Raton,  Clayton,  and  Capulin.  The 
average  free  air  anomaly  over  the  map  is  about  +55  mgal,  with  a  large 
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free  air  anomaly  of  +70  mgal  associated  with  Capulin  Mountain.  This 
large  free  air  anomaly  at  the  mountain  station  (+70  mgal)  is  +15  mgal 
different  than  the  regional  free  air  value  immediately  around  the  moun¬ 
tain  (+55  mgal).  The  335  m  cone  comprising  the  mountain  has  approximately 
the  same  volume  as  a  right  cylinder  112  m  high  with  the  same  radius  of 
500  m.  The  gravitational  effect  of  such  a  cylinder,  assuming  a  density 
of  2.67  g/cc,  would  be  12.4  mgal  if  the  top  of  the  cylinder  is  located 
at  mountain  peak  elevation.  This  calculated  value  is  not  precisely  equal 
to  the  15  mgal  anomaly  observed  associated  with  Capulin  Mountain,  but 
the  method  of  calculation  is  only  approximate.  Agreement  is  close 
enough  to  allow  the  assumption  that  the  gravitational  effect  of  Capulin 
Mountain  itself  is  15  mgal. 

Bouguer  Gravity  Data 

The  plot  of  average  elevation  against  average  Bouguer  anomalies 
(from  figure  6)  indicates  an  inverse  relationship  between  elevation  and 
the  Bouguer  data  (figure  9) .  The  Bouguer  anomaly  is  actually  the  nega¬ 
tive  gravity  effect  of  the  root  increment  below  the  compensation  level. 

As  the  terrain  elevation  increases,  the  root  increment  below  the  compen¬ 
sation  level  also  increases  in  thickness  to  maintain  isostatic  equilibri¬ 
um.  An  increase  in  the  thickness  of  the  root  below  compensation  level 
results  in  a  more  negative  Bouguer  anomaly.  The  Bouguer  gravity  maps 
(in  hand  drawn  form  -  figure  10,  and  computer  plotted  form  -  figure  11) 
Indicate  this  relationship  between  rising  topography  (figure  7)  and 
a  decreasing  Bouguer  equipotential  surface.  The  slab  assumed  in 
making  the  Bouguer  correction  results  in  the  computation  of  a  Bouguer 
anomaly  that  is  excessively  negative  in  areas  which  are  distinguished 
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Figure  9.  Plot  of  the  averaged  NOAA  elevation  data  vs,  the  averaged 

NOAA  Bouguer  gravity  data.  The  slope  of  the  line,  m,  shows 
the  inverse  relationship  of  elevation  and  Bouguer  gravity 
values . 
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Figure  10.  Regional  Bouguer  gravity  map  using  NOAA  gravity  data. 

The  dots  represent  the  gravity  stations.  Datum  is 
sea  level,  and  a  density  of  2.67  g/cc  was  used  in  cal¬ 
culating  the  Bouguer  correction. 
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Figure  11.  Regional  Bouguer  anomaly  map  of  the  study  area.  Gravity  sta¬ 
tions  are  the  same  as  in  figure  10.  Contour  interval  is  10 
mgal,  datum  is  sea  level. 
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by  topographically  isolated  terrain,  such  as  mountains.  Sierra  Grande 
(fountain  is  an  example  of  this  effect;  however,  Spanish  Peaks  and  the 
Vet  (fountains  regions  have  sufficient  areal  extent  to  allow  the  slab 
assumption  in  the  Bouguer  calculation. 

Eagles  Nest  (fountain  has  associated  with  it  both  a  positive  free  air 
anomaly  and  positive  relief  with  respect  to  the  surrounding  Bouguer  sur¬ 
face  on  the  Bouguer  anomaly  map.  The  positive  free  air  anomaly  is  ex¬ 
pected  due  to  the  near  surface  mass  effect  of  the  mountain.  The  presence 
of  positive  relief  on  the  Bouguer  equipotential  gravity  surface  may  indi¬ 
cate  a  near  surface  high  density  body,  which  would  cause  the  Bouguer  cor¬ 
rection  to  be  too  small,  or  it  may  indicate  the  extension  below  datum  of 
the  igneous  material  forming  the  mountain.  In  either  case,  it  appears 
that  the  density  of  the  rocks  forming  Eagles  Nest  Mountain  is  greater 
than  2.67  g/cc. 

The  free  air  anomalies  corrected  for  the  Faye  effect  in  the  study 
area  have  already  been  shown  to  indicate  isostatic  adjustment,  and  the 
computation  of  isostatic  anomalies  directly  from  Bouguer  anomalies  can  be 
used  to  add  validity  to  this  conclusion.  The  computation  of  the  isostat¬ 
ic  anomalies  done  here  is  based  upon  the  Airy  isostatic  model  with  the 
density  contrasts  and  formulas  suggested  by  Woollard  (1969) .  The  slab 
thickness  of  the  mass  below  the  compensation  level  is  determined  by 
assuming  that  part  of  the  negative  Bouguer  anomaly  is  caused  by  the  pres¬ 
ence  of  a  slab  of  crustal  material  displacing  heavier  mantle  material. 
Woollard  (1969)  determined  that  a  crustal  density  of  2.92  g/cc  is  repre¬ 
sentative  of  the  crustal  rock  at  the  compensation  level  (M-discontinuity) 
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and  a  density  of  2.67  g/cc  is  representative  of  the  upper  crustal  rock. 
The  gravitational  effect  of  the  elevated  surface  is  equal  to  the  simple 
Bouguer  correction  computed  using  a  density  of  2.67.  However,  the  effect 
of  the  surface  slab  upon  a  measured  value  on  the  surface  is  exaggerated 
as  compared  to  the  effect  of  the  compensating  root  slab  upon  the  same 
measurement  as  previously  discussed.  Therefore,  if  the  region  is  in  iso¬ 
static  balance,  the  observed  simple  Bouguer  anomaly  is  equal  to  the 
Bouguer  gravity  effect  computed  for  the  surface  slab  with  average  eleva¬ 
tion  in  the  area* minus  the  exaggeration  effect  of  slab  position.  Pre¬ 
viously,  it  was  shown  that  this  exaggeration  effect  is  approximately 
equal  to  the  observed  free  air  anomaly.  The  isostatic  anomalies  for  each 
one  degree  block  in  the  study  area  have  been  calculated  in  this  manner, 
and  the  values  are  shown  in  table  1. 

These  isostatic  anomalies  show  the  study  area  to  be  in  near  isostat¬ 
ic  equilibrium,  as  was  suggested  in  the  free  air  anomaly  discussion. 
Qureshy  (1962)  studied  the  Bouguer  data  for  all  of  Colorado  and  concluded 
that  the  eastern  portion  of  the  state,  including  the  northern  part  of 
this  report's  study  area,  is  in  isostatic  equilibrium.  The  slight  posi¬ 
tive  bias  in  the  results  may  be  explained  by  either  crustal  thickening 
or  increased  crustal  density  (Woollard,  1969). 

Table  1  includes  calculations  of  crustal  thickness  using  Shurbet's 
(1966)  density  model  and  formulas.  The  crust  thins  from  west  to  east, 
and  continues  to  thin  east  of  the  study  area.  Stewart  and  Pakiser  (1962) 
calculated  a  seismic  refraction  profile  using  the  Gnome  explosion. 

Their  computed  thickness  is  very  close  to  the  crustal  thicknesses 


calculated  using  gravity  data  in  this  report. 

In  addition  to  the  NOAA  Bouguer  gravity  data,  a  separate  Capulin 
Bouguer  gravity  map  (figure  12)  was  prepared  from  data  collected  by  the 
author.  The  Bouguer  anomaly  associated  with  Capulin  Mountain  is  about 
-200  mgal,  which  is  approximately  20  mgal  more  negative  than  the  adjacent 
area.  Since  the  mountain  is  a  feature  with  little  area,  the  slab  assumed 
in  the  Bouguer  correction  has  exaggerated  the  Bouguer  anomaly.  Because 
the  mountain  itself  is  335  m  high,  a  total  of  37.4  mgal  was  subtracted 
from  the  measured  value  as  a  part  of  the  total  Bouguer  correction  for  the 
height  of  the  mountain.  The  actual  gravitational  effect  of  the  mountain 
has  previously  been  shown  to  be  about  15  mgal.  Therefore,  the  Bouguer 
correction  applied  to  the  mountain  station  was  exaggerated  by  22.4  mgal. 
Correction  for  this  exaggeration  shows  the  area  beneath  the  mountain  to 
be  only  slightly  more  negative  than  the  surrounding  area,  and  no  ad¬ 
ditional  mass  deficiency  exists  beneath  the  mountain  as  compared  to  the 
surrounding  area. 

Since  the  study  area,  including  the  Capulin  area,  has  been  shown  to 
be  in  isostatic  balance,  the  Bouguer  equipotential  surface  is  expected  to 
reflect  variations  in  crustal  thickness  as  well  as  local  changes  in 
crustal  density.  Therefore,  the  Bouguer  anomalies  can  be  used  to  cal¬ 
culate  sedimentary  thickness  in  the  Basin.  Calculations  of  this  kind 
were  carried  out  to  produce  a  basement  relief  or  structure  map. 

Raton  Basin  Basement  Structure  Map 

In  previous  studies  (Griggs,  1948;  Baltz,  1965),  subsurface 
structure  in  the  Basin  has  been  deduced  from  surface  rocks  and  well  log 


Figure  12.  Capulin  Bouguer  gravity  map.  Data  stations  are  the  same 
as  shown  in  Figure  8.  Density  used  for  the  Bouguer  cor¬ 
rections  is  2.67  g/cc. 
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data.  A  basement  relief  map  based  on  structural  data  cannot  be  drawn 
because  geologic  mapping  in  the  Basin  is  incomplete.  However,  Bouguer 
gravity  data  were  used  in  this  study  to  deduce  crustal  thicknesses  and 
to  determine  basement  structure  (or  relief) .  The  basement  surface 
computed  is  defined  as  the  boundary  between  the  less  dense  sedimentary 
rocks  above  and  denser  crystalline  rocks  below.  Intrusive  rocks  that 
penetrate  the  strata  are  also  classified  as  basement  rocks  because  their 
presence  causes  effective  thinning  of  the  sedimentary  column  in  the 
computation. 

The  best  method  for  determining  basement  relief  is  based  upon  the 
realization  that  the  Bouguer  anomaly  is  actually  the  negative  gravity 
effect  of  the  slab  or  root  which  has  displaced  the  heavier  material  of 
the  sub-compensation  level.  In  this  study,  it  is  assumed  that  the  neg¬ 
ative  effect  is  locally  affected  by  less  dense  sedimentary  rocks  dis¬ 
placing  heavier  basement  rocks;  therefore,  the  variation  in  the  Bouguer 
gravity  values  can  be  interpreted  in  terms  of  variation  in  the  thickness 
of  the  sedimentary  rocks. 

As  a  first  step  in  constructing  a  basement  relief  map  for  the  Basin 
it  was  necessary  to  ascertain  the  regional  trend  of  the  Bouguer  gravity 
map.  A  second  order  polynomial  surface  was  empirically  determined  to 
best  represent  the  Bouguer  gravity  surface  as  discussed  in  Chapter  II. 
Variation  between  this  regional  trend  and  the  observed  local  Bouguer 
gravity  values  is  assumed  to  be  caused  by  differences  in  thickness  and 
density  of  the  sedimentary  rocks. 

The  areal  extent  of  the  various  lithologic  units  is  shown  in  plate 
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1;  Clark  (1966)  gives  representative  densities  for  the  rock  types  pres¬ 
ent.  Representative  densities  for  the  columnar  sections  (figure  3)  for 
the  Las  Vegas  and  Raton  Basins,  respectively,  were  found  to  be  2.27  g/cc 
based  on  the  densities  reported  by  Clark  (1966) .  The  density  contrast 
between  the  representative  density  for  the  Basin  and  the  density  used  in 
the  Bouguer  correction  (2.67  g/cc)  is  -0.4  g/cc.  Estimations  of  lateral 
variations  in  the  density  contrast  throughout  the  Basin  were  based  on 
the  following  criteria:  1)  sedimentary  rock  thinning  near  the  Apishipa, 
Cimarron,  and  Sierra  Grande  Arches;  2)  presence  of  volcanic  rock  in  the 
southern  portion  of  the  Basin;  3)  higher  density  intrusive  rocks  near 
Spanish  Peaks  and  Wet  Mountains  regions;  and  4)  an  average  sedimentary 
rock  density  of  1.8  g/cc  (producing  a  -0.87  g/cc  density  contrast)  for 
the  Great  Plains  (east  of  the  Apishipa  Arch  -  Sierra  Grande  Arch  axes) 
as  determined  by  Shurbet  (1966) .  A  representative  calculation  from  the 
Bouguer  gravity  map  (figure  10)  along  37°N,  104°-105°W  is  included 
(figure  13) . 

At  37°N,  104. 5°W  (figure  13),  the  Bouguer  anomaly  value,  determined 
from  the  Bouguer  gravity  map,  is  -216  mgal.  The  regional  Bouguer  trend 
value  (-210  mgal)  was  computed  from  a  second  order  polynomial  surface 
described  in  Chapter  II.  The  residual  gravity  value  (Ag)  at  37°N, 

104. 5°W  is  the  difference  between  the  Bouguer  gravity  value  and  regional 
trend  value  at  that  point  (-210  +  216  mgal  ■  +6  mgal) .  The  density  con¬ 
trast  value  (Ap)  at  37°N,  104. 5°W,  determined  from  the  density  contrast 
model,  is  -0.58  g/cc.  The  depth  to  the  basement  relief  surface  is  found 


Figure  13.  Basement  structure  profile  derived  from  NOAA  Bouguer  grav¬ 
ity  data  along  37°N,  104°-105°W. 

A.  The  regional  Bouguer  dip  is  determined  by  a  smooth 
curve  through  the  Bouguer  anomaly  equipotential  surface. 

B.  The  difference  of  the  regional  Bouguer  dip  from  the 
Bouguer  anomaly  equipotential  surface  in  A  is  known  as  the 
residual  gravity  (Ag) . 

C.  The  density  contrast  model  (Ap)  is  constructed  from 
the  known  densities  of  the  stratigraphic  units  in  the 
Basin  and  subtracting  those  densities  from  the  2.67  g/cc 
density  used  in  the  Bouguer  correction.  Since  the  sedi¬ 
ment  density  is  1.8  g/cc  and  it  is  assumed  that  the  sedi¬ 
ments  are  displaced  crustal  rock  (p  ■  2.67  g/cc),  the 
gravity  effect  of  this  displacement  is  negative. 

D.  The  depth  to  the  basement  structure  surface  (h)  is 
calculated  by  using  Ag  from  B  and  Ap  from  C  above  in  the 
gravity  mass  effect  formula:  h  -  Ag/2irApY. 
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by  solving  the  gravitational  effect  formula  for  depth  (h),  using  the  Ag 
and  Ap  values  above:  h  »  Ag/2iTYAp  ■  -810  feet  (below  sea  level).  This 
type  of  calculation  was  repeated  for  each  of  the  2,400  density  contrasts 
in  the  Basin,  and  the  resulting  depths  to  basement  structure  were  then 
contoured  by  the  computer  techniques  described  in  Chapter  II.  The  ac¬ 
curacy  of  the  depths  to  the  basement  structure  surface  was  independently 
verified  using  well  log  data  (Baltz,  1965)  and  a  measured  section  from 
the  Raton  Basin  (McGehee,  1955). 

In  order  to  describe  the  correlation  between  the  geology  on  plate  1 
and  the  basement  structure  on  plate  2,  the  same  topographic  names  are 
used  to  refer  to  features  on  both  maps.  This  practice  is  justified  be¬ 
cause  the  features  named  on  the  basement  structure  map  all  have  assocated 
topographic  features.  In  the  area  of  the  Wet  Mountains  and  Spanish 
Peaks,  the  basement  structure  map  shows  that  basement  rock  is  exposed; 
that  is,  there  are  no  sedimentary  rocks  present.  The  structurally  high 
Apishlpa  Arch  ana  on  the  basement  relief  map  extends  southeast  of  the 
Wet  Mountains.  The  basement  rises  at  about  1.5°  to  approximately  1,500 
feet  of  vertical  relief  over  a  ten  mile  distance.  The  Apishlpa  Arch 
area  is  outlined  by  positive  relief  for  approximately  50  miles  on  the 
basement  structure  map.  The  basement  structure  map  shows  the  Apishlpa 
Arch  as  the  eastern  boundary  to  the  Basin. 

The  Delcarbon  and  La  Veta  syncline  areas  are  shown  on  the  basement 
structure  map  as  depressions.  The  Delcarbon  depression  is  above  sea 
level  and  the  La  Veta  depression  is  below  sea  level;  both  depressions 
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deepen  to  the  south.  The  axis  of  the  La  Veta  depression  is  the  deepest 
depression  on  the  basement  structure  map  and  defines  the  Raton  Basin 
axis. 

A  positive  relief  feature  on  the  basement  relief  map  between  the 
La  Veta  and  Delcarbon  depressions  is  also  associated  with  intrusive  rocks 
on  the  geologic  structure  map  (plate  1) .  The  shallow  depth  to  known  oil 
bearing  formations  of  Cretaceous  age  (Baltz,  1965)  make  this  anticline  a 
favorable  site  for  petroleum  exploration.  However,  the  heat  from  the 
intrusive  activity  may  have  destroyed  any  oil  present,  in  a  manner  simi¬ 
lar  to  the  destruction  of  coal  beds  found  elsewhere  in  the  Basin  and  des¬ 
cribed  in  Chapter  I. 

The  Cimarron  Arch  is  an  area  of  positive  basement  relief  and  clearly 
separates  the  Raton  and  Las  Vegas  Basins.  It  is  asymmetric  in  shape, 
steeper  on  the  north  slope  (A. 5°)  than  the  south  slope  (2.9°),  and  has 
approximately  3,500  feet  of  vertical  relief.  The  Cimarron  Arch  area  ex¬ 
tends  about  fifty  miles,  from  the  west  edge  of  the  basement  structure  map 
to  the  Sierra  Grande  Arch  area.  The  marked  relief  of  this  basement  fea¬ 
ture  suggests  the  Arch  has  effectively  prevented  intertonguing  of  Las 
Vegas  and  Raton  Basin's  sedimentary  rocks.  As  the  Basins  were  downwarped 
and  isostatic  adjustment  occurred  to  compensate  for  the  accumulating  sed¬ 
iments,  the  Cimarron  Arch  remained  a  positive  structural  feature. 

Although  not  a  significant  topographic  or  gravity  feature  today,  the 
Arch  was  important  in  the  development  of  the  Basin. 

The  Sierra  Grande  Arch  area  is  roughly  defined  by  moderate  relief 
on  the  basement  structure  map  for  approximately  twenty  miles  east  of 
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the  Las  Vegas  Basin.  Total  positive  relief  increases  about  1,000 
feet  over  a  twelve  mile  distance,  producing  a  gentle  0.8°  slope.  The 
whole  eastern  sector  of  the  basement  structure  map  is  a  low  relief 
positive  area,  rising  slowly  to  the  Apishipa  and  Sierra  Grande  Arch 
areas . 

A  previously  unmapped  and  unnamed  basement  trough  which  is  in  near 
alignment  with  a  left  lateral  fault  south  of  Spanish  Peaks  (plate  1), 
trends  northeast  from  Raton  Basin  toward  Apishipa  Arch,  The  trough  deep¬ 
ens  southeast,  towards  the  Raton  Basin,  and  may  be  fault  controlled.  As 
a  structurally  low  feature,  the  unnamed  trough  may  contain  significantly 
more  sediments  than  areas  surrounding  it.  The  margins  of  the  unnamed 
trough  may  prove  favorable  for  future  petroleum  exploration.  Presence 
of  this  trough  is  confirmed  on  the  Raton  Basin  magnetic  map. 

Raton  Basin  Magnetic  Map 

The  Raton  Basin  vertical  magnetic  data  is  tabulated  and  included  in 
Appendix  C.  The  magnetic  data  is  shown  in  both  hand  drawn  (figure  14) 
and  computer  plotted  formats  (plate  3) .  The  computer  plotted  magnetic 
map  is  drawn  at  the  same  scale  as  the  basement  structure  map  for  com¬ 
parison.  Zeitz  et  al.,  (1969)  showed  that  the  Colorado  sedimentary 
rocks  in  the  Basin  are  essentially  nonmagnetic,  the  metamorphic  and  vol¬ 
canic  rocks  are  weakly  to  moderately  magnetic,  and  the  Precambrian 
quartz  monzonite  rocks  are  moderately  to  strongly  magnetic.  The  Raton 
Basin  magnetic  maps  actually  represent  the  surface  of  the  basement  rocks 
in  the  Basin.  Magnetic  anomalies  may  also  be  associated  with  lateral 
changes  in  rock  types,  but  there  is  no  geologic  evidence  for  basement 
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rock  type  changes  in  the  Basin  with  the  exception  of  those  areas  of 
known  intrusive  rocks. 

Metasedimentary  rocks  are  exposed  near  the  Cimarron  Arch  and  these 
rocks  are  near  station  number  one,  which  was  the  magnetic  survey  datum. 
Magnetic  measurements  in  the  area  of  the  metasedimentary  rocks  gave  no 
indication  that  these  rocks  have  higher  magnetic  susceptibilities  than 
the  sedimentary  rocks  in  the  Basin.  Therefore,  it  is  reasonable  to 
assume  that  the  rocks  in  New  Mexico  have  the  same  magnetic  susceptibil¬ 
ities  that  Zeitz  e£  al^. ,  (1969)  reported  for  Colorado  rock  types.  Some 
lack  of  detail  in  the  magnetic  map  in  the  southern  portion  of  Raton 
Basin  is  the  result  of  poor  road  access.  All  magnetic  measurements  were 
made  along  roads. 

Structural  depressions  (plate  1)  are  associated  with  negative  mag¬ 
netic  anomalies  (plate  3)  in  the  Basin.  The  La  Veta  and  Delcarbon  syn¬ 
clines  both  have  distinct  low  magnetic  anomalies  associated  with  them; 
the  unnamed  trough  detected  on  the  basement  structure  map  also  has  a  dis¬ 
tinct  low  magnetic  anomaly  associated  with  it.  The  magnetic  anomaly  as¬ 
sociated  with  the  unnamed  trough  makes  it  appear  to  deepen  to  the  north¬ 
east,  opposite  the  direction  determined  on  the  basement  structure  map. 
However,  this  apparent  deepening  probably  reflects  only  the  computer 
response  to  the  lack  of  magnetic  stations  in  the  eastern  part  of  the  sur¬ 
vey  area. 

Positive  magnetic  anomalies  are  associated  with  positive  structural 
features.  For  example,  the  unnamed  anticline  between  the  La  Veta  and 
Delcarbon  synclines  has  a  significant  positive  magnetic  anomaly 


associated  with  it.  The  Cretaceous  age  sediments  with  known  oil  shows 
are  only  thinly  buried  near  this  anticline. 

Capulin  vertical  magnetic  data  were  also  studied  as  an  addition  to 
the  Raton  Basin  magnetic  data.  The  Capulin  data  are  tabulated  and  in¬ 
cluded  in  Appendix  B  with  the  Capulin  gravity  data.  Unlike  the  sediment 
filled  Raton  Basin,  the  Capulin  area  is  covered  by  volcanic  rocks,  which 
have  significantly  higher  magnetic  susceptibilities  than  the  sedimentary 
rocks  beneath.  Some  of  the  magnetic  anomalies  in  the  Cauplin  area 
(figure  15)  are  associated  with  topographic  features.  For  example,  there 
is  a  magnetic  anomaly  associated  with  Capulin  Mountain.  However,  a  posi¬ 
tive  magnetic  anomaly  on  the  east  slope  of  the  mountain  is  not  associated 
with  a  topographic  feature.  Free  air  gravity  anomalies  also  mark  this 
area,  and  this  anomaly  relationship  may  be  an  indication  of  lateral  vents 
or  dikes. 

Variation  in  the  magnetic  properties  of  different  volcanic  rock 
types  (Dobrin,  1976)  would  be  helpful  in  identifying  magnetic  anomaly 
patterns  in  volcanic  regions.  Indeed,  the  largest  positive  magnetic 
anomaly  is  near  Sierra  Grande  Mountain;  smaller  but  equally  distinctive 
positive  magnetic  anomalies  are  near  Capulin  Mountain.  The  unusually 
high  magnetic  susceptibility  of  the  Sierra  Grande  sequence  (indicated  by 
its  magnetic  anomaly)  may  be  diagnostic  of  these  volcanic  rocks  and  aid 
in  determining  their  areal  extent. 

Capulin  Local  Seismic  Study 

In  addition  to  gravity  and  magnetic  measurements,  an  effort  was 


made  to  determine  natural  seismicity  in  the  Capulin  area.  It  was 


Figure  15.  Capulin  magnetic  map.  The  dots  represent  magnetic 

stations.  Station  number  one  is  used  as  reference  and 
is  represented  by  a  star. 
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thought  that  some  tectonic  activity  may  be  in  progress.  The  re¬ 
sults  of  the  seismic  survey  in  the  Capulin  area  and  all  known  seismic 
activity  during  the  Twentieth  Century  are  included  in  table  2,  The  epi¬ 
centers  of  the  earthquakes  in  table  2  are  also  shown  on  a  map  of  north¬ 
eastern  New  Mexico  (figure  16) . 

A  bar  graph  (figure  17)  indicates  the  distribution  of  seismic 
activity  in  northeast  New  Mexico  from  1900  to  1980.  The  lack  of  any 
significant  recorded  seismic  activity  during  the  first  half  of  this 
century  is  due  largely  to  the  lack  of  recording  equipment  or  populace 
to  observe  any  shocks.  Since  1962  the  number  of  recorded  events  have 
increased  due  to  the  installation  of  permanent  seismograph  networks. 

In  the  last  five  years  seismic  activity  in  New  Mexico  has  declined 
(Sanford,  1981),  yet  this  study  suggests  that  the  Capulin  area  has 
rather  high  seismic  activity.  The  events  recorded  during  this  study  were 
local  in  origin  and  small  in  magnitude,  and  were  not  recorded  elsewhere 


Figure  16.  Geographic  location  of  seismic  events  strong  enough  to  be 
felt  by  local  population  in  northeast  New  Mexico  from 
1900  to  1980. 
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Table  2.  Recorded  Earthquake  Activity  in  the  Study  area  of  New  Mexico 
1900-1981. 


SOURCE 

DATE 

YR-MO-DAY 

TIME 

GMT 

LAT 

°N 

LONG 

°W 

INTENSITY  modified  Mercali 
or  local  magnitude 

1 

06-04-19 

2100 

V 

1 

06-04-20 

1300 

- 

- 

VI 

2 

24-08-13 

0423 

36.0 

104.5 

V 

2 

52-08-13 

2042 

36.5 

105.0 

V 

2 

63-06-06 

0806 

36.6 

104.4 

2.7 

2 

66-09-24 

0734 

36,4 

105.1 

2.7 

2 

66-09-24 

0827 

36.4 

105.1 

2.4 

2 

66-09-25 

1011 

36.3 

105.1 

2.8 

2 

66-09-25 

1223 

36.5 

105.1 

2.8 

2 

72-02-20 

2310 

36.4 

104.9 

1.5 

2 

72-02-20 

2323 

36.4 

104.9 

2.2 

2 

75-05-16 

0138 

36.9 

105.0 

1.5 

2 

75-05-16 

0726 

36.5 

104.7 

1.9 

2 

75-06-21 

0542 

36.1 

104.0 

2.0 

3 

79-09-19 

0539 

36.4 

103.7 

1.8 

3 

80-09-12 

2138 

36.5 

105.1 

2.3 

4 

81-07-24 

1613 

4 

81-07-24 

1815 

4 

81-08-27 

0220 

4 

81-08-27 

0238 

4 

81-07-27 

1554 

4 

81-07-27 

1738 

4 

81-07-30 

1824 

4 

81-08-25 

2034 

4 

81-08-27 

2041 

4 

81-09-01 

2124 

4 

81-09-01 

2254 

4 

81-09-09 

1724 

4 

81-09-16 

0336 

Source  Data: 

1.  Anonymous,  The  Raton  Range  (April  21,  1906)  shock  felt  at  Folsom,  NM 

2.  Sanf ord  ejt  al . ,  (1981) 

3.  Sanford  (1981) 

4.  Capulin  data  (this  report)  epicenters  assumed  at  or  near  Capulln,  NM. 
No  intensities  were  calculated,  but  the  sensitivity  of  the  equipment 
would  indicate  the  relative  intensities  to  be  much  smaller  than  pre¬ 
viously  recorded  in  the  area. 


1950  1960  1970  1980 

EARS) 

f  seismic  events  In  Table  2  in  five  year 
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in  the  state  (J.J.  Wolff,  personal  communication,  1981).  Some  of  the 
events  recorded  may  have  been  local  mining  blasts  or  sonic  booms, 
although  sonic  booms  are  usually  easily  distinguished  (figure  18) . 

The  presence  of  seismicity  near  Capulin  suggests  this  area  to  be 
either  affected  by  dormant  vulcanism,  or  by  tectonic  adjustment  along 
faults  in  the  area.  Sanford  at  al . ,  (1981)  suggest  that  seismic  activity 
in  northeast  Mew  Mexico  may  be  due  to  an  eastern  extension  of  the  Jemez 
Lineament  in  western  New  Mexico.  Continued  monitoring  of  the  Capulin 
seismicity  would  be  necessary  to  determine  the  source  and  extent  of  the 
local  seismic  activity. 


CHAPTER  IV 


CONCLUSIONS 

The  geophysical  measurements  and  the  interpretation  reported  in  this 
thesis  meet  the  objectives  set  out  in  Chapter  I.  Furthermore,  the  in¬ 
terpretation  of  the  data  of  these  geophysical  surveys  lead  to  the  follow¬ 
ing  results: 

1.  The  free  air  and  Bouguer  gravity  data  both  indicate  the  study 
area  is  in  near  perfect  isostatic  adjustment.  Correlation  between  this 
report  and  studies  in  adjacent  areas  Indicate  the  Colorado  -  northeast 
New  Mexico  -  west  Texas  region  is  in  isostatic  balance. 

2.  A  large  free  air  gravity  anomaly  in  the  southern  portion  of  the 
study  area  is  the  result  of  local  crustal  loading  by  Quaternary  volcanic 
activity.  The  thickness  of  the  volcanic  rock  is  estimated  to  180  m 
based  on  the  gravitational  effect  represented  in  the  free  air  anomaly. 

3.  The  basement  relief  map  constructed  from  Bouguer  gravity  data 
has  a  high  correlation  with  known  geologic  structures. 

4.  The  Raton  magnetic  map  generally  represents  the  same  basement 
surface  as  that  constructed  from  Bouguer  gravity  data. 

5.  The  Capulin  area  and  adjacent  High  Plains  are  seismically 
active.  Further  investigation  of  the  local  seismicity  is  important. 

6.  The  basement  structure  map  and  Raton  magnetic  map  indicate  areas 
of  potential  benefit  in  future  petroleum  exploration.  The  high  heat  flow 
in  the  Basin  may  have  economic  importance.  If  the  heat  flow  is  found  to 
be  a  result  of  dormant  vulcanism,  seismic  studies,  as  in  the  Capulin 

area,  may  be  helpful  in  determining  the  extent  of  the  thermal  field. 
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APPENDIX  A 

NOAA  GRAVITY  AND  ELEVATION  DATA 


Explanation  of  Table  Headings 

Lat  Northern  latitude  of  station  in  degrees. 

Long  Western  longitude  of  station  in  degrees. 

Elev  Station  elevation  in  meters  above  sea  level. 

FAA  Free  air  anomaly,  calculated  from  a  datum  of  sea  level.  The 

free  air  correction  (FAC)  is  0.3085  mgal/m  and  is  applied  to 
the  observed  gravity  reading  (OGR) :  FAA  *  OGR  +  FAC. 

BA  Bouguer  anomaly,  calculated  using  a  datum  of  sea  level,  station 

elevation,  and  a  density  of  2.67  g/cc.  The  Bouguer  correction 
(BC)  is  27ryph,  or  0.0837  mgal/m,  and  is  applied  to  the  FAA: 

BA  -  FAA  -  BC. 


L-T  J--’ 


■U-Ei. 
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1  _  4.  1  -.3 

3 a.  3 1  o 

104. 3oG 

bo.  6  1  9 

1  a  4. Goh 

3  a  .  5  3  e 

i ’4.29  t 

JO.  4  J  V 

1-4.  -*E  J 

3  o.  6 1  ft 

lot.  30  5 

jo.  /  E  2 

1  -  *7  .  1  J  O 

ft  a.  7  a  ft 

1  ft  -t .  48  a 

i  o  .  ft  4  2 

1-4.  jCft 

3o.519 

10 4. e54 

EO.  7  O'  3 

1  0  4.  7 <7-  2 

fto. 39  4 

1  o  4.  44  3 

3  o. ai 9 

1-  ■*♦.  a  _  e 

3a. 713 

lC  4.  5<+C 

bo.  C.24 

L  -  4  .  ft 1  ( 

3  0.53-4 

lU  <7.  5*.— 

?  6. ft  5  1 

1  -  4  •  Oft  -E 

36. 593 

104.  737 

3  0.  ft4  E 

L.,’4.  219 

36.  o  1  2 

104.915 

<o,57? 

1  2  4  .  9  a  2 

3a.3Ea- 

1 J4.  5t  7 

J  —  .  ft  O  J 

10  4.7  3 ft 

ft  5  .  o  2  0 

1 C  t .  b  7  ft 

jo.  o  e  3 

1-4.573 

3  o. so  7 

1C  4.  79  2 

3o.7  37 

1  -  4  •  ft  9  G 

3  a.  .>0  8 

1C 4. 71 2 

3  a.  J  5 2 

10  4.  5  4  f 

3a.  39  7 

134.755 

j  a  .  ft  6  o 

1-  4.  9ft  ft 

ft  6 .  *5  6  9 
30. 5  J4 

3ft. 536 
3  c  .  o'*  3 
3  a.  a  9  C 
bo.  i'  I  h 
3ft.  75  ti- 
3  o  .  o  t  o 
3ft.  ^  j  •— 

3 o.  33j 
3a. 51 0 

jo.  bi! 

2  a.  cs't 
"3  o .  *•  ~  ? 

3  ft. 

ft  o . 

3ft. 

*>  • 

3a. 
bo. 

3a. 

ft  o  • 

bo. 
ft  a. 
ft  a. 

bo. 


Eft  2 
aac 
osi 
^  ft  ? 
it  9 
a39 
9  /  9 
■J  J4 

*1  3 
921 
4  _e 
a2ft 
7ft  9 


lift. 79o 
1  840 

1ft  *♦.  7fto 
104. 37J 
i :  *♦.  fttE. 
lCt.9C4 
1  -  4  •  9  2  1 
l  ift-.  8ft  7 
l-ftft.  337 
lo4.  909 
i  -.1  4.  ft  ft  j 
lCft. 793 
134. ?). 
lift.  7  75 

i C  4 . toft 

lift. 62  5 
l  ...b  ft  2 
1  Oft. 3b9 
l  ft. 937 

10ft. 971 

1 3  4  •  9  e  5 

104.376 

1  -  ft.  7c  ft 

lift.VOft 

1  v.'  ft.  )ft2 


.  ft 

■  V 


1  ft  ?  a  < 

1  936 , 

1  700  I 

2  1 2 1  •  t  -  C 
2 31 o .ova 

»_E.Eo.aoG 
?229.ftCG 
2  1 14.700 
22.-j.i00 
1955.0CG 

EOiiEa  /  E 

2  lb  6 . kcc 
2  2-2.  a  -  G 
Zlfte.ft'.'O 

.JE  1  •  ftG-U 

2  tt  ft  •  bo  G 
2279.460 
2136 . 3CC 
a  e29.  e ^  iy 
25*o. 220 
e  e o .  GOO 
2ft  3  3  .  fty  G 

tftll.oi-a 

2179. CCG 

2  'n’.jaJ 


21io»-*cc 

5  c  .  t  G  C 

2  115.  -  -  v> 

ft  A  .  ft  C  E 

£ x o* .320 

ft  a  .  4  -  C 

2‘2-l.ftCO 

E  7 . 0  C  e 

2C8G. auC 

39.7C0 

—  O  EC  .  —  -  0 

e  ft . 9  3 E 

2  C  0  a  .  4C  0 

33. ICC 

9  I4  7 . 3  C  - 

4b.  V  -  c 

21C8.2- J 

46. 81 G 

-  0  E4.5ou 

44  •  A  2-  » 

2COC.uoO 

2ft.  bCC 

2  1e>6. 5v. 0 

4*2  .  A  1  6 

l9ftl.ftCC 

1 3  .  ft  C  a 

l? aft. 2 CO 

0 1. 4O  e 

19  78.-.  '  J 

2 1 . 8  i.  e 

---3.-12 

A  ■  .  7C  C 

1  7S1?.  -E  A 

22. ?CU 

2  1 1  7  .  1  ■  0 

ft  /  .  Y  c  E 

1421. eCO 

.  f*  L  - 

5  -J4. b-G 

ft  •«  •  4 2  2 

22-l.aoO 

fta  .  -CC 

_0j4.  ftA  0 

21.  i  a'  G 

2C34.5oo 

2ft.  ICC 

1  /_  1  .  JoC 

,3,oi  - 

1 9  C  2 . 0  0  0 

lo.2C  - 

1 0— 5 .oi¬ 

i  •  ft  C  G 

ly-  7  .ft  CO 

7  *  -  — 

1  -  -  3  .  e  —  a 

—1  »  ft  c  c 

1677. eoO 

1 5  .  ft  C  G 

1  94-V  .3'  u 

1 0 .  /  e  c 

1  cab . aOC 

12. fta  C 

1 40- .  f c  0 

.  EEC 

1625. 2'-0 

5.52  v’ 

1  4  e  4  •  0  —  0 

a2 . -G- 

1  c?9 3 . 72u 

Aft . 460 

I  >34.  ?'-u 

3  .  ft  E  y 

1 9ftC . 5CC 

1  a  .  <5  e  G 

1444. 4E- 

aa.lv- 

2 199.500 

27.0-C 

E  1  4  4 . 4C  G 

i  j.’<e  e 

2429 . 4-00 

.a  i.  COG 

2aE,o.E-a 

?.vvC 

2C8C • bCC 

12.  iCC 

-i to . -V  o 

-184.-  >/l 
-1 aft.2Co 
— loo . ftCG 

—19  ft.CGC 
-iGo.vCi. 


A  4  »  E  OE 

11.8vG 
19.  IV  0 

26. COG 
37. ft ft o 

31. 9Go 
_-i.  ICC 
1ft.  1  >  -  o 

AO. 9G i 


.1  o 


ICC 
-19b. 2C0 
-lBb.fttC 
-ltft.vU 
-193 .  3CC 
— ifto. foe 
— _02 • oCC 
-19ft. iOC 
-lftft.ftUL 
— 207 ioOe 
— 14  ft .oC  0 
-lftv..Cl 
-206. -CC 

—  1 6  C .  ECO 
-19 C. ICC 
-20  5.  t Co 
— _ 0  *♦  .  ftC  c 

— l9l . A-'C 

—  1 9  6 • ftoC 
-19  7 . 1  Co 
-195.  ICC 
— 20  i  .  ->Oo 
-194. jGG 
-2 01. ICe 

—206 .7U4 

-20v. ftCC 
-196 .aCC 
— eC  e  .  3o  v. 

—  .9  ft . ftOO 
-1 9  -ft  . 1 C C 
-199. 3Go 
— 19ft . -Co 
— 21 o .ftCO 
—2  2 1 . ft  0  C 
-233. 3CC 
-214.  e ft  C 
—2 1 4 . 9  C  C 
— 2C2.  e>CC 
—2  Go .oi4 
-200.  .OC 
-211. -0C 
-221. ftCO 
—21  •  ftCC 
— 2 1  a  •  j  C  C 
-217.  bCC 
-23a;.  3Cc 


10 . -C  0 

—20  3 . hCG 

-a  z.  zee 

14. 7  e  e 

— e29. 7CC 

■>  •  ft  G  r. 

—  2  3  9 . 9  C  v 

9.aC0 

-242. zee 

i  9  .  j  (_'  L 

-64 1 .-CC 

37. ICC 

—240  .-CO 

9  .  eC  0 

* ■  w  4  j  •  c  U/ 

jO.tCv 

-206  ,-C  - 

2  7 . oC C 

-iftft. ?Cv 

3?  .2'1C 

-<:5o.-Cv 

it . oe  e 

-249. jCC 

3O.2C0 

-E47. ICC 

ft4  •  L'tO 

— e  4  — . 32  - 

11 .OCC 

-2ft  ft. lee 

1  -'C  - ‘ 

-224.  - CC 

»m**>  t  *  ■-  ■m 


I 

u 


i  f 

L0  6L 

..  .  =  4A 

5  A 

3t>»  aC  5 

1 0  8£C 

2166. 5C0 

19.700 

—221 .4-00 

3  :>.  C  '  7 

1  -•  4  •  4  I  1 

2562.500 

2  7 . 3  0  w 

-236.500 

36.029 

U4.ooa 

1421 ,4LG 

21.340 

-143.500 

55.5  12 

10  4.77*. 

1  ‘#00 . 40  0 

-i.-CC 

-612. >C0 

3o. 62  b 

1 o4. 756 

19C7 .100 

19 .100 

-194.100 

3c. o2  5 

1 .4.523 

1  9fc  7.300 

t3.4CC 

-190.3CC 

30.61*t 

104. 943 

2146. 40  C 

26.4C0 

-215.700 

50. 054 

1 0  4.  9fc  7 

2  523,  joo 

38 .eOc 

-221. ICC 

36. 5o7 

106.651 

1970. 20  0 

15.3GC 

—205. wCO 

3  o.  v>3  7 

i04. o3o 

1799. 600 

3  .  **C  0 

-197.700 

3o. 666 

104.633 

1 3o 7 .200 

7.3C0 

-201.500 

36.31  i 

104.  92  7 

1  9  6  4 . 5  c  0 

10.100 

— <1 1  O’  •  1  w  0 

56.  316 

10  4.  946 

1-59  L.o(  0 

4.9CC 

-212. 7oC 

50. 5i 6 

104.  9e*+ 

2cofc  «  act' 

9 .200 

-215.500 

3o.52  3 

10  4.497 

2039 .4t  0 

0. 4OO 

-221. OCC 

3  5.  3o  7 

104. 926 

23 . uCC 

-2x0 .4CO 

3  o. 117 

104. 660 

1916. uGo 

I4.0OO 

-20l.oOC 

3o. 1 1 7 

1  6  65 1 

19C j.uCO 

7.60  0 

—40'  5  •  nO C 

50.117 

10  4.  56  6 

1  90  2  •  ‘70 G 

7.400 

— 2  0  5 . 5  0  C 

3o. 115 

12  6.423 

2 19  i.ico 

35.800 

-204. _cC 

36.114 

10  4.961 

2254. uCO 

31.2C0 

-220. 9GC 

3  o. 115 

i. 0 4.  5i2 

16  lb.  /CO 

0  .iv.  O' 

— L9o.6cC 

3o. 11 3 

10-..  771 

bOti. 40 0 

25.600 

-202.500 

36. 3o2 

10  6. 593 

1 723.200 

4.  CCl 

— 164.0 Ou 

3  6.00  7 

106.  7CC 

1 b69. 3CC 

1.10C 

— 210.20*0 

3  o  .  1 V  3 

1-4.  660 

1023.40. 

4  •  4  0  4 

-144.200 

36.  2t?  6 

1  3  4. od  5 

1909 .  -0  0 

4.  cCc 

-209.400 

36.095 

1  v 4-1 .  7o  7 

L  tC  4  .  30  J 

4 . 7 ;  c 

-2Co.5CC 

36.  113 

1C4.S90 

2150.300 

35 .500 

-205. OCC 

3  O.  lot 

1  j  4  .  59  1 

I4O5.40O 

12.500 

-20O.7oc 

36.177 

10  4.  y9b 

2 169. oCO 

19.4-00 

-223.3CC 

50.235 

10**.  6  55 

1955.600 

12.106 

—  2  o  .  4  C  C 

3o.^2§ 

1  v  4.  69  5 

1830 . oOO 

4.9CC' 

-200. 7CC 

5  6  .  2  i  5 

i'>,  6y  0 

1  64  5  .  -»0  0 

3.5*0 

-2(iC.“v  0 

50.  3  l  J 

*  j  **•  3*5 

2  -*.S  •  oc  0 

51.700 

-197. ^CC 

3  o. 32  5 

lv*..  92  0 

1.I6. 1.0 

1  6 , 4  O  C 

--Co  .  uv-’O 

o  6, 36  3 

104.717 

162  >. 5  00 

5.400 

— 19  4  • GO 

3  o.  4  J  6 

1  w  4  .  6  1  7 

ib39.20o 

4  •  54  0 

-201. 2C0 

3q.41 3 

136. a5« 

1915.900 

7. 30  0 

— 20a. oC 0 

3c.. 477 

104.75/ 

1  j 7 4. 60  O 

-7.5C0 

-21  7.500 

36.026 

106.5  37 

1137. 3C0 

23.500 

-137.600 

3  6.066 

1  .  -»  •  o9  5 

I.eO.oou 

I.4CO 

—208 . 4uC 

36. 39  7 

106.937 

2C50. 7C0 

22.100 

-207. 2CC 

3o.  «.oS 

1 0  4. 8 66 

1923. 0C0 

3.55c 

-204.500 

36.46C 

106.950 

2018. ICO 

lo. ICC 

-209. 5CC 

36. 33  0 

104.923 

1990. 500 

1 4 . 6  b  C 

—20  7.4(  0 

3  6. 3o  o 

104.397 

1940. 10 0 

12.200 

-204. ICO 

36.294 

10  6.  92  7 

2071 . 1  CO  • 

2 0 . 30  c 

-211.500 

3o.  117 

1C  4. 666 

1913.500 

6  .  *t0  0 

-205. 700 

56. 117 

1^/4.  5V6 

1864..  6C0 

1 . 500 

-20 4.5C0 

36.4  55 

IC6.843 

1902.300 

C  •  Oo  C 

-204.100 

50.456 

104.792 

<.  5a2  •  0C0 

v.2i'C 

-210. 2C0 

3o. 480 

I04.77o 

1373.900 

-6. ICO 

-215.700 

3  O  .  i  J  5 

1  -  4  .  b.4  jf 

2_Oo. -Cu 

46. 700 

-20c. 700 

56. 264 

12  4.  o37 

1416.100 

13. 40  C 

-200.  71  0 

3  o  •  2  o  1 

1  J  4  .  7  6  5 

1354.400 

1 0  .  7  6  0 

-20C.1C0 

36.  346 

106.792 

1836 .500 

15.500 

— l95.oC0 

5  6.0  j  6 

1 0  4 . 6  5  C 

1447.100 

i:.40C 

-20  7.200 

Jo  .0  76 

104.730 

1932.400 

2.400 

-207.100 

50.071 

1 0  4.  99  3 

2 1 40 »  uOo 

25 . 700 

-219.10c 

36.066 

104.592 

1917.200 

10.600 

-203.3C0 

50. 056 

104. 52/ 

1 870. 60  J 

17 .300' 

-191. 8CC 

36.044 

1  w  4, 6  32 

1928.503 

9  »a0  0 

—  206. wCO 

3o. 0  5  C 

1C  4. 73 x 

200-1 .  *-00 

i.  7  .  cOG 

—  2  0  6  *  0  C  0 

56. 02 0 

lu4. 3lo 

1995.300 

13. COO 

-2O5.4OO 

3  6. o2  2 

1.4. 7 5t 

l'/ti5.;L0 

11.300 

-20b.2cc 

5  6 . 0  1  5 

1C  4,  546 

1 WC  1 «  3cC 

<-3.400 

-l£S .ocC 

3  o.  Ol  5 

1'.  4.  a5u 

1426.500 

8 . 5uC 

-20  i.xCC 

Copy  available  to  DTIC  does  sol 
penrnt  fully  l.gible  tii production 
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3o.  0 1  3 
3  a.  Ci  2  j 
3o.023 
3o. C  2  7 
36.034 

3  6, 0  ti 

36.  3  i  0 
3o.  1  ■+  3 
3o. CDS 
3o.  •*-»o 
2o.477 
3o.-*°4 
36. 29u 
j>6.  2“  C 
3o. 27  6 
3o. 136 
36.197 
3  o  .  0  a. 

5o.  15  9 
3  b.  17  3 
36. IS  1 
3 o.  5“t 
36.  33  C 
36.39  5 
36. 0u6 
36.Ho 
3  o  .  1  6 
3o. „3  5 
3o.3o2 
3b.3  7  3 
3o.  3  3  <_ 
36.317 
3o.273 
36.  23“ 
3o. Ia5 
36. 1 j  7 
36. 100 
5  6.3b! 

30.44  3 

3  6.  -*-9  « 

36. 45  $ 
3  o .  1 3  0 
36.42.£ 
3  6  .  “  1 9 
3  o  .  4  J  C 
3o. 3v  7 
3o.377 
3o. 330 
3o.i50 
3o. 3bC 
3o. 2 1 o 
3o. 31o 
3b. 29C 
3  b, 44  b 

3o»  23  o 
36.  2  7 “ 
3  o  .  2  3  0 
io. 21 7 
36.216 
3o. 1 3 2 
3o.  179 
36.  l-»0 
3b. 0  b  2 

3o.“23 
3o.323 
3  o  .  2  b  3 
3o.  24  0 
5  o«  1  “  3 


1)4.614 
1 j“.44“ 

1  J  4.  7“1 

i  v'u 
104.  8oo 
10-*.  924 
l  J3.  it  5 
Lb2. Ibi 
10J.2C3 
1  „  6  .  “4  o 
lc  3. 260 
10  3.  3b 2 
1  33.  “t>3 
1  w  3 . 40  b 
Iu3.  27o 
10  3. “76 
103.439 
10 3. 3b“ 
103.41  / 
lC 2.275 
1 0  3»  3““ 
1.  3. 455 
1 0  3  .  2b  3 
10 3. 3“ 5 
1 0  3  •  4g  o 
1 Z 3. “26 
K 3. 49 S 
102. 3k3 
103.433 

10  5*“bu 

10  3.029 
10  2 . C“  1 
103.041 
103. 941 
193.040 
10  3.041 
103.040 
103.02/ 
133. 023 
103. 233 
103.020 
103.177 
103. 152 
103.020 
103. 066 
i. 3  .  1 6  v 
LC  3 • 0^6 
103.  173 
103.125 
1C3.C70 
103.  160 
102. 013 
103.  132 
10  3.07  i 
1)5. ISO 
1  b  3  .  w  “  b 
10  3.03a 
103. 125 
103.062 
103.03c 
1 j5. 132 

lO  3.05/ 
103.037 
1C  3. “50 
103.310 
10  3. “96 
1  3  •  J4  6 

103.127 


194-..  00  0 

£  4,  **  G  •  W  G  O 

1933 . 200 
2 15b.  2'- 0 
2090. oOO 
2  1  3  5  .  b‘5  b 
1476.690 
1  “  i  o  .  7Co 
1407. 500 
a /2  9-  j 

1 535.900 
1  to  1 .7  CO 
1603. “00 
ItCC.oCO 
1 3l“. 9b  0 

1 5  6 1 . 2  o  0 
1535 .000 

1 “6 1 . btC 
1521.610 
l4c5.bOC 
L4o7 . OCO 

loli.o.-J 
1650. 300 
1  57  ^  .  30’j 
1 4c  3 .600 
1 4o2 . ICO 
i“7“. 000 
1550.3-0 
1637. ICO 
1 5  'C  4  •  40  0 
1414. oOO 
1  *tC  b  .  o0  C 
1346.500 
1  3  n  j  .  4 ..  3 
1374. jOG 
1351 .200 
13o / .2u0 
l“2i * “b0 
1439.  --GG 
1 “2 5. 4-0 

14^.6 . 800 
1 53o.600 
1  50  7.  I i/G 
l  >04.  7CC 
1459. 100 
l  “9  5  .  0 

1 4<_4  .  lC  0 
147b •  bOv 
1459 .“00 
1  “3  7  .  /  C  0 
l“7c.  ICO 
i399.„-0 
l“4o • oCO 
i “ l 3 • “00 
1“ 7 7.200 
l“CC.“Cw 
13c2.300 
l“33.6w0 
1412. LOO 
1  >3“.  “OC 
14C5. iGO 
1 3o / • bCO 
13ol»200 
1 b“b . “C  3 
1 625.700 
io“l»7c0 
142 1  *  “i’0 
1 j9c. IuC 


13.400 
_  /  .  bt  C 

9. >cc 

37.400 
32. CC0 

50.400 

“1  .lbU 
-ll.lbb 

5  «  9CG 
4  7 .  t>C  C 

12.  CG0 

-  J»V  b  b 

13 . 3b C 
1  6  »  “C  0 

1 . 3b0 
1 4  .  o  5  0 

13.  “C  0 
3.200 

11 . 700 

1 .  o  0  C 

4.  200 
17. *CC 

6. 10C 
3  .  3C  b 
2o.  c.C  0 
li . bt 5 
1  o  .  7  0  C 
“1.7.0 
15.100 
7.70b 
-1.300 
-1. ICO 
7CC 
—4  .“t  v 
—3.70  0 
—2 . 00  0 
-9. 300 
“  .  C  b  C. 

5. “0C 
•  o  0  0 

l.oOO 
7. 10C. 

3  •  bCC 
—0  .—CO 
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Station 

Lat 

Long 

Free  Air  An. 

Bouguer  An. 


CAPULIN  GRAVITY  AND  MAGNETIC  DATA 


Explanation  of  Table  Headings 

Numerical  listing  of  gravity  and  magnetic  station  iden¬ 
tification.  Stations  with  apostrophe  (')  are  addition¬ 
al  magnetic  stations. 

Northern  latitude  of  station  in  degrees. 

Western  longitude  of  station  in  degrees. 

Free  air  anomaly  (FAA)  calculated  from  a  datum  of  sea 
level.  The  free  air  correction  (FAC)  is  0.3085  mgal/m 
and  is  applied  to  the  observed  gravity  reading  (OGR) : 
FAA  -  OGR  +  FAC. 

Bouguer  anomaly  calculated  using  a  datum  of  sea  level, 
station  elevation,  and  a  density  of  2.67  g/cc.  The 
Bouguer  correction  (BC)  Is  2Tryph,  or  0.0837  mgal/m 
and  is  applied  to  the  FAA:  BA  *  FAA  -  BC. 

Relative  vertical  magnetic  reading  in  gammas.  Datum  is 
station  one. 


Rel.  Mag. 


1 

36.805 

103.956 

49.0 

-186,2 

0 

2 

36.806 

103.951 

48.8 

-183.1 

695 

2' 

36.804 

103.960 

1710 

3 

36.800 

103.965 

57.89 

-184.5 

1395 

3' 

36.796 

103.963 

2070 

4 

36.792 

103.963 

64.0 

-179.5 

925 

4' 

36.790 

103.965 

850 

5 

36.794 

103.961 

59.25 

-185.8 

1340 

5' 

36.792 

103.966 

1790 

6 

36.790 

103.978 

60.0 

-187.1 

1950 

7 

36.796 

103.982 

58.5 

-186.6 

275 

7’ 

36.798 

103.979 

-1200 

8 

36.809 

103.980 

53.0 

-183.5 

750 

9 

36.785 

103.983 

57.85 

-189.3 

75 

10 

36.782 

103.969 

57.06 

-188.7 

145 

11 

36.776 

103.973 

54.4 

-187.9 

600 

12 

36.767 

103.969 

50.6 

-186.7 

275 

12* 

36.755 

103.978 

-150 

13 

36.759 

103.979 

49.3 

-186.6 

200 

14 

36.752 

103.978 

54.8 

-181.1 

450 

15 

36.743 

103.978 

47.63 

-186.3 

1125 

16 

36.743 

103.992 

47.8 

-186.2 

1500 

17 

36.744 

103.984 

47.7 

-186.2 

2425 

17’ 

36.735 

103.979 

1250 

18 

36.744 

103.973 

47.9 

-186.1 

2125 

19 

36.746 

103.962 

47.5 

-184.5 

2350 

20 

36.746 

103.951 

49.05 

-182.6 

2775 

21 

36.748 

103.938 

49.1 

-181.6 

3625 

22 

36.748 

103.927 

50.05 

-181.2 

3525 

23 

36.749 

103.917 

49.9 

-181.3 

3500 

24 

36.779 

103.979 

58.45 

-189.4 

1375 

25 

36.784 

103.964 

59.3 

-204.2 

3500 

26 

36. 783 

103.971 

74.6 

-195.0 

3800 

27 

36.756 

103.917 

49.85 

-180.2 

4650 

28 

36.759 

103.927 

48.7 

-180.6 

4600 

28' 

36.762 

103.926 

4525 

29 

36.764 

103.927 

49.1 

-180.5 

4500 

30 

36.764 

103.917 

49.5 

-179.7 

4750 

30' 

36. /61 

103.917 

5000 

31 

36.769 

103.917 

48.1 

-179.3 

4700 

32 

36.742 

103.936 

49.9 

-184.2 

4920 

33 

36.734 

103.936 

52.5 

-182.5 

3750 

34 

36.727 

103.935 

52.3 

-183.2 

6050 

35 

36.723 

103.935 

51.6 

-184.5 

5500 

36 

36.716 

103.935 

51.3 

-185.5 

3650 

37 

36.713 

103.935 

50.85 

-185.7 

3860 

38 

36.713 

103.944 

50.3 

-184.7 

3950 

39 

36.713 

103.950 

48.8 

-185.0 

46  .0 

72 


Station 

Lat 

.  Lcmg _ 

Free  Air  An. 

Bouguer  An. 

Rel,  Max 

40 

36.705 

103.953 

48.0 

-185.4 

3300 

41 

36.698 

103.952 

47,75 

-185.0 

3600 

42 

36.698 

103.945 

47.8 

-185,4 

4000 

43 

36.698 

103.934 

49.3 

-185.1 

3600 

44 

36.741 

103.968 

49.1 

-184.5 

3450 

45 

36.741 

103.958 

48.9 

-183.3 

3600 

46 

36.732 

103.958 

48.0 

-184.5 

2450 

47 

36.727 

103.953 

48.0 

-184.4 

3750 

48 

36.719 

103.953 

48.4 

-184.7 

3500 

49 

36.689 

103.965 

48.4 

-185.3 

2850 

50 

36.689 

103.972 

48.6 

-185.6 

3200 

51 

36.689 

103.988 

48.2 

-186.4 

2700 

52 

36.689 

104.023 

48.1 

-187.7 

2000 

53 

36.689 

103.999 

48.3 

-186.5 

2350 

54 

36.712 

103.998 

52.1 

-186.3 

1750 

55 

36.721 

103.998 

49.3 

-187.7 

2000 

56 

36.727 

103.999 

48.5 

-186.9 

1650 

57 

36.734 

103.998 

47.9 

-188.6 

1900 

58 

36.739 

103.996 

47.8 

-187.3 

1900 

60 

36.743 

104.003 

47.3 

-186.1 

1375 

61 

36.743 

104.011 

45.9 

-186.9 

1700 

62 

36.743 

104.021 

45.85 

-187.7 

1275 

63 

36.746 

104.042 

44.85 

-188.3 

1375 

64 

36.727 

104.007 

47.7 

-186.7 

1225 

65 

36.726 

104.016 

46.3 

-187.5 

1400 

66 

36.726 

104.025 

47.7 

-187.5 

1725 

67 

36.708 

104.025 

45.7 

-188.9 

1425 

68 

36.708 

104.036 

46.8 

-189.8 

1325 

69 

36.787 

104.005 

52.2 

-188.2 

1775 

70 

36.784 

104.000 

52.55 

-187.6 

1375 

71 

36.781 

104.001 

51.4 

-188.1 

1500 

72 

36.828 

103.990 

53.9 

-189.8 

1725 

73 

36.789 

103.958 

54.4 

-187.2 

5075 

74 

36.777 

103.962 

54.4 

-186.0 

4000 

75 

36.777 

103.974 

57.9 

-188.3 

2675 

76 

36.792 

103.966 

54.4 

-185.7 

2125 

77 

36.786 

103.960 

68.1 

-180.6 

3900 

78 

36.789 

103.964 

58.1 

-188.8 

2975 

79 

36.815 

103.945 

49.5 

-178.4 

2250 

73 
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MAG 


APPENDIX  C 

RATON  BASIN  MAGNETIC  DATA 

Explanation  of  Table  Headings 
Northern  latitude  of  station  in  degrees. 

Western  longitude  of  station  in  degrees. 

Relative  vertical  magnetic  reading  in  gammas.  Station  one  is 
datum. 
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PLATE  I 

GEOLOGY  AND  STRUCTURE  OF  THE 
RATON  BASIN  OF  N.M.  AND  COLO. 
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PLATE  2 

BASK M KMT  STRUCTURE  MAP  KOR  TUK 
RATON  HAS! N  OK  N.M.  AND  COLO. 
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PLATE  3 
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